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Abstract — In  this  paper  a method  for  the  solu- 
tion of  scattering  problems  is  proposed.  In  partic- 
ular the  EFIE  is  solved  by  a collocation  point  pro- 
cedure, where  the  unknown  current  density  is  ex- 
panded by  Daubechies  wavelets  on  the  interval  and 
the  integration  is  performed  in  the  wavelet  domain 
by  the  integral  operator,  hence  without  the  use  of 
any  quadrature  formula.  Comparison  with  induced 
currents  calculated  by  a standard  MoM  and  with 
fields  calculated  by  a FEM  code  are  reported. 


I.  Introduction 

Wavelet  Expansion  (WE)  has  become  a widely 
used  tool  in  electromagnetic  analysis.  The  main 
reasons  can  be  found  in  the  strong  interpolating 
properties  of  the  wavelet  functions  and  to  the  fact 
that  WE  and  reconstruction  can  be  performed  by 
the  use  of  fast  algorithms  (see  for  example  [l]  - [3]). 

The  analysis  of  scattering  problems  can  be  car- 
ried out  expanding  the  unknown  functions  in  terms 
of  a chosen  wavelet  basis  and  performing  a Galcrkin 
procedure  using  the  same  wavelet  functions  as  test 
functions  (sec  [4]  - [6]).  The  integration  is  then 
performed  by  the  use  of  quadrature  formulae. 

In  this  paper  the  EFIE  for  the  scattering  on  a 
conductive  body  is  solved  in  a different  way.  By  uti- 
lizing Daubechies  wavelets  on  the  interval  the  un- 
known current  is  expanded,  then  a collocation  point 
method  is  used  and  the  integration  is  performed 
by  the  use  of  the  integral  operator  for  Daubechies 
wavelets  on  the  interval  developed  by  the  authors 
(see  [8]).  In  this  way  the  need  of  quadrature  formu- 
lae is  avoided;  furthermore  the  well  known  capabil- 
ity of  the  wavelet  functions  of  representing  irregular 
signals  with  few  coefficients,  allows  us  to  use  bases 
of  low  dimension  (if  compared  with  the  number  of 
unknowns  of  a standard  MoM);  for  these  reasons 
both  accuracy  and  CPU  time  saving  are  achieved. 


II.  Wavelets  on  the  Interval  and 
Operators 

The  concepts  of  scaling  functions,  wavelets,  time- 
scale  analysis,  multiresolution  analysis  are  here 
considered  known  [1];  there  are  many  wavelet 
bases  available  in  the  literature,  and  we  chose  the 
Daubechies  Wavelets  on  the  interval  [2]  for  their 
numerical  properties.  In  particular  the  choice  of 
wavelets  that  ’’survive”  only  on  intervals  is  adopted 
because  we  arc  interested  in  the  solution  of  a bound- 
ary value  problem. 

From  the  wavelet  theory  we  know  that  WE  must 
be  performed  starting  from  a signal  known  at  a 
dyadic  number  of  samples.  This  number  is  equal  to 
the  dimension  of  the  basis  on  which  we  perform  the 
WE;  hence  a signal  represented  by  n = 2 771  sam- 
ples when  expanded  in  the  wavelet  domain  leads 
to  a number  of  2m  coefficients.  For  the  com- 
pact support  wavelet  there  is  an  important  rela- 
tion that  allows  a straight  computation  of  the  co- 
efficients at  the  higher  resolution  of  a generic  func- 
tion: it  is  possible  to  obtain  them  from  the  sam- 
ples of  the  functions  itself  according  to  the  relation 
(<t>j,k,f)  = 2J/2f(2Jk),  where  <pj,k  is  the  scaling 
function  of  order  J,  k of  the  adopted  wavelet  basis. 
Then,  the  vector  of  wavelet  coefficients  g represent- 
ing the  wavelet  transform  of  a function  g(x)  can  be 
obtained  by  multiplying  a matrix  W related  to  the 
adopted  wavelet  basis  and  the  samples  g(xj)  corre- 
sponding to  2m  equally  spaced  points  in  the  inter- 
val [0,1].  Further  details  about  wavelet  numerical 
computation  can  be  found  in  [3]. 

When  performing  WE  of  a function,  the  nota- 
tion that  will  be  used  throughout  the  paper  is  the 
following: 

/(<)  = = b(a:)f  (!) 

i 

where  b(x)  = [6* (a;),  • • • 6n(x)]  is  the  wavelet  basis 
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and  f = [/i,--*,/n]T  is  the  vector  of  coefficients 
constituting  the  wavelet  expansion  of  the  signal. 

Also  operators  can  be  represented  in  the  wavelet 
domain,  as  described  in  [7].  The  authors  have  ob- 
tained the  representation  of  the  integral  operator 
for  the  Daubechies  wavelets  on  the  interval  (see  for 
example  [8]).  The  convenience  in  using  operators  in 
the  wavelet  domain  is  that  the  integral  of  a function 
/(x)  can  be  calculated  by  the  matrix-vector  prod- 
uct Tf  where  T is  the  constant  sparse  matrix  rep- 
resenting the  operator  in  the  wavelet  domain  and 
f is  the  wavelet  expansion  of  the  function  f(x).  In 
the  previous  operation  the  result  is  the  primitive 
function  of  /.  In  order  to  better  understand  the 
previous  statement  let  us  define  f(x)  = x3;  we  can 
write  that 

f f(x)dx  = f x3dx  — ~x4  = g(x)  (2) 
[0,1]  [0,1] 

where  the  integral  limit  [0, 1]  indicates  that  we  are 
considering  the  interval  [0,1].  Expanding  (2)  ac- 
cording to  (1)  in  the  wavelet  domain  we  can  write 
that 

f b(x)f dx  — b(x)g  (3) 

[o,i] 

where  f and  g are  the  wavelet  expansion  of  the  two 
functions.  Left  multiplying  by  b(x)T  and  taking 
into  account  the  definition  of  the  integral  operator 

( T = (b(x)  ,/ b(x)))  and  the  orthogonality  prop- 
erties of  the  wavelet  basis  (la  = (b(x),  b(x)),  whith 
la  being  the  identity  matrix)  we  can  write  that 

Tf  = g (4) 

Hence  as  clearly  stated  in  equation  (4)  the  cal- 
culation of  function  g(x)  can  be  performed  in  the 
wavelet  domain  by  multiplying  matrix  T by  vector 
f,  the  wavelet  expansion  of  the  function  /(x);  then 
the  result,  vector  g,  is  inverse  transformed,  obtain- 
ing the  function  g(x). 

In  order  to  compute  a definite  integration  directly 
in  the  wavelet  domain,  two  new  (row)  vectors  must 
be  introduced:  vectors  bo  = b(0)  and  bx  = b(l) 
that  are  the  values  of  the  basis  functions  on  the  left 
and  right  border,  respectively,  of  the  interval  [0, 1]. 
The  meaning  of  these  two  vectors  is  the  following: 
given  a function  /(x),  its  left  border  value  /( 0)  can 
be  obtained  by  the  coefficients  of  its  WE  simply  by 
the  use  of  the  vector  b0  as  /(0)  = b0f. 


Suppose  that  we  want  to  compute 

l l 

J f(x)dx  = J x3dx  = ^x4\l  = p(a:)|i  = i (5) 

0 0 

Then  starting  from  equation  (4)  it  must  be  con- 
sidered that  the  primitive  function  (g)  is  already 
calculated,  and  we  only  need  to  evaluate  it  on  the 
borders  of  the  interval.  Hence  it  yields 

l 

J f(x)dx  = big  - b0g  (6) 

0 

which  can  be  rewritten  as 

1 

J f(x)dx  = b,Tf  — b0Tf  (7) 

0 

As  evidenced  in  equation  (7)  the  calculation  of 
the  definite  integral  of  a function  /(x)  can  be  per- 
formed knowing  its  WE  f , the  integral  operator  ma- 
trix T and  the  border  vectors.  The  quantities  T, 
bo  and  bi  are  known  once  a wavelet  basis  has  been 
chosen,  so  they  need  to  be  computed  only  once,  and 
not  at  any  analysis. 

III.  Method  of  Moments  and  Wavelet 
Expansion 

A.  General  Considerations 

In  the  study  of  scattering  from  conducting  cylin- 
ders, an  integral  equation  can  be  formulated,  which 
in  general  has  the  form  of 

\p2  J j Js(r')G(rs,r')ds'+ 

+V  //  V'  • Js(r')G(rs,  r')ds'  = E‘(r  = rs) 

S (8) 

where  7/  = y/(fi/e)  and  (5 2 = Js  is  the  cur- 

rent density  induced  on  the  scatterer,  G is  the  green 
function  for  the  three  dimensional  scatterer,  r'  and 
rs  are  respectively  the  integration  variable  and  the 
observation  point,  both  on  the  surface  scatterer; 
and  EJ  is  the  incident  field. 

In  a simpler  way  and  in  one  dimension  equation 
(8)  can  be  in  general  rewritten  as 

J g(xf)K(xixf)dxf  + c{x)g{x)  = h(x)  (9) 

where  g(x)  is  our  unknown  function.  In  the  lit- 
erature [5],  function  g is  expanded  in  the  wavelet 
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domain,  and  the  expansion  is  substituted  in  (9).  In 
order  to  obtain  a linear  system  for  the  unknown 
coefficients,  the  resultant  equation  is  tested  with 
the  same  expansion  functions  (Galerkin’s  method). 
Quadrature  formulae  available  in  the  literature  for 
wavelets  and  scaling  functions  are  used,  and  a 
square  system  is  obtained,  which  solved  gives  the 
unknown  coefficients. 

Due  to  the  high  interpolating  properties  of  the 
Daubechics  wavelets  on  the  interval,  we  have 
adopted  this  family  of  wavelets,  for  which  the  au- 
thors have  developed  the  representation  of  the  inte- 
gral operator.  Wavelets  on  the  interval  (and  so  op- 
erators in  the  wavelet  domain)  are  defined  on  [0, 1], 
hence  the  contour  of  the  scattercr  must  be  mapped 
into  the  interval  [0,1].  As  it  is  suggested  in  [6]  for 
an  arbitrary  contour  of  the  scattcrer  two  steps  must 
be  performed: 

• The  contour  of  the  scattcrer  is  discretized  in 
boundary  elements  and  then  each  boundary  ele- 
ment (simply  a first  order  element)  is  mapped  into 
one  dimensional  standard  element  through  shape 
functions  or  interpolation  functions. 

• The  standard  elements  are  mapped  into  corre- 
sponding portions  of  the  interval  [0,1]. 

In  this  way  the  basis  functions  arc  defined  in  a 
standard  way  on  the  interval  [0, 1],  since  the  contour 
has  been  mapped  on  this  interval. 

B.  Scattering  from  a Conductive  Body 

In  case  of  a two  dimensional  problem  with  TMZ 
polarization  the  EFIE  equation  is  the  following: 

Y f Up')H^(l3\pm  - p'\)dc'  = Ei(pm)  (10) 

where  pm  is  any  observation  point  on  the  scattercr, 
p'  is  any  source  point  on  the  scattcrer  and  C is  the 
contour  of  the  scattercr. 

After  the  mapping  on  the  elemental  interval  [0, 1] 
(described  in  the  previous  section)  is  performed,  we 
perform  a classical  collocation  point  procedure:  we 
evaluate  equation  (10)  at  a particular  point  pm  of 
the  contour;  hence  equation  (10)  can  be  rewritten 
as 

l 

k f }{x)g{x)dx  - h (11) 

o 

where  k = f(x ) = H^\f5\pm  - p'|)  and  g(x)  - 

Jz(p'). 

Then  the  wavelet  expansion  in  the  space  domain 
is  performed.  It  has  to  be  noticed  that  in  equation 


(11)  there  is  the  product  between  two  functions,  in 
particular  one  of  them  is  known  ( f(x ))  while  the 
other  is  our  unknown.  As  explained  in  [8]  it  is 
possible  to  obtain  the  wavelet  expansion  of  a the 
product  between  two  functions  y(x)  = f(x)g(x)  as 
the  product  between  a constant  diagonal  matrix  F 
and  a vector  g,  where  g is  the  wavelet  expansion  of 
the  function  g(x)  and  F is  a diagonal  matrix  whose 
entries  arc  the  samples  at  n = 2m  equally  spaced 
points  of  the  values  of  f(x)  in  the  interval.  This 
approximation  is  as  much  accurate  as  the  number 
of  samples  is  high,  hence  as  the  resolution  of  the 
chosen  basis  increases.  In  case  the  two  functions 
are  know  the  above  described  procedure  is  useless, 
since  the  expansion  of  y(x)  can  be  performed  di- 
rectly. But  in  the  cases  when  one  of  the  two  func- 
tions is  unknown  then  the  procedure  is  fundamen- 
tal, since  it  allows  to  keep  g as  the  unknown  vector 
and  anyway  perform  the  wavelet  expansion.  Hence 
equation  (11)  can  be  expressed  in  the  wavelet  do- 
main as 

l 

k J b(x)Fgdx  = b(a;)h  (12) 

o 

where  F is  the  diagonal  matrix  with  the  samples 
of  the  Hankel  function,  g is  the  vector  of  unknown 
coefficients  and  h is  the  expansion  of  constant  h 
on  the  interval  [0, 1]  i.e.  considered  as  a constant 
function  on  the  whole  interval.  By  left  multiplying 
equation  (12)  by  b(x)T  we  obtain 

l 

kb(x)T  J b(x)dxFg  = b(x)Tb(x)h  (13) 
o 

and  taking  into  account  the  definition  of  the  inte- 
gral operator  and  the  orthogonality  properties  of 
the  chosen  basis  we  obtain 

/cTFg  - h (14) 

In  this  way  equation  (14)  establishes  a relation 
between  the  primitive  of  the  product  f(x)g(x)  and 
the  constant  function  h(x)  = constant.  This  is 
actually  something  different  from  what  we  want, 
hence  the  two  vectors  bo  and  bi  (introduced  in  sec- 
tion II)  must  be  employed.  Hence  we  can  write  that 

bifcTFg  - bo/cTFg  = h (15) 

Equation  (15)  is  characterized  by  known  matri- 
ces T and  F and  known  vectors  boand  biand  by 
the  unknown  vector  of  coefficients  g.  Once  the  res- 
olution of  the  wavelet  basis  is  chosen  the  number  of 
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basis  functions  is  consequently  fixed,  hence  equa- 
tion (15)  must  be  written  for  a number  of  n points 
on  the  interval  itself.  This  leads  to  a sparse  linear 
square  system  whose  unknowns  are  the  coefficients 
g and  which  can  be  solved  in  low  CPU  time. 

IV.  Numerical  Results 

The  numerical  results  presented  here  are  relative 
to  the  scattering  of  a square  conductive  object  illu- 
minated by  a polarized  TMZ  field.  The  geometry 
of  the  system  and  the  input  signal  are  reported  in 
figures  1 and  2. 

a y 


Fig.  1.  Geometry  of  the  system 


calculated  current  density  on  the  surface  scatterer 
for  a number  of  64,  128  and  256  wavelet  functions 
for  a frequency  of  1 GHz. 


Fig.  3.  Current  density  on  the  scatterer  calculated  by  the 
use  of  64  wavelets 

It  is  evident  that  together  with  the  increase  of  the 
resolution  from  64  to  256,  the  accuracy  of  the  re- 
sults becomingh  higher.  Nevertheless  even  at  lower 
resolutions  the  obtained  current  (at  a very  lo  CPU 
time  cost)  is  consistent  witht  the  problem. 


Fig.  4.  Current  density  on  the  scatterer  calculated  by  the 
use  of  128  wavelets 


Fig.  2.  Input  signal 

The  diagonal  of  the  scatterer  is  of  0.2m  and  the 
frequency  content  of  the  input  signal  is  of  the  order 
of  GHz.  FFT  has  been  used  in  order  to  obtain  the 
behavior  in  the  time  domain.  Figures  3-5  show  the 


Figure  6 shows  the  comparison  between  the  cal- 
culations performed  by  the  proposed  method  with 
a resolution  of  256  wavelets  and  a standard  MoM 
(with  collocation  point)  technique  with  500  points 
on  the  whole  perimeter.  It  can  be  seen  the  very 
good  agreement  between  the  two  different  methods. 
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Fig.  5.  Current  density  on  the  scatterer  calculated  by  the 
use  of  256 


Fig.  6.  Current  density  on  the  scatterer  calculated  by  the 
use  wavelets  and  MoM 

As  for  the  CPU  time  the  method  proposed  is 
of  the  order  of  2-4  times  faster  than  the  standard 
MoM,  due  to  the  high  numerical  efficiency  of  the 
wavelet  expansion  and  integration,  as  described  in 
the  previous  sections:  the  construction  of  the  inte- 
gral operator  matrix  and  the  border  vectors  is  done 
only  once  and  can  be  seen  as  a pre  processing  activ- 
ity, while  the  integration  is  performed  by  a simple 
matrix  - vector  product,  without  the  need  of  any 
quadrature  formula.  Furthermore  the  well  known 
numerical  properties  of  the  wavelet  functions  (well 
addressed  in  the  literature)  allow  the  choice  of  bases 
of  small  dimensions  in  order  to  obtain  accurate  re- 
sults. 

Since  the  calculation  of  the  current  density  is  not 


the  only  important  result  from  an  engineering  poin 
of  view,  a comparison  with  results  obtained  by  the 
use  of  a FEM  code  of  the  total  field  are  also  re- 
ported. In  figures  7 and  8 the  electric  field  in  point 
P = (-0.2,0)  and  P = (0,  — 0.2)  are  reported, 
evaluated  by  the  FEM  code  and  by  the  proposed 
method,  by  using  a wavelet  basis  of  256  functions. 


Fig.  7.  Electric  field  evaluated  at  point  P = (—0.2,0) 


Fig.  8.  Electric  field  evaluated  at  point  P = (0,  —0.2) 


Figure  9 reports  the  comparison  between  two  so- 
lutions obtained  by  the  proposed  method  at  dif- 
ferent resolutions,  in  particular  with  64  and  256 
wavelets,  and  shows  the  robustness  of  the  method 
in  terms  of  calculated  fields 
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Fig.  9.  Comparison  at  different  resolutions  of  Electric  field 
evaluated  at  point  P = (0,  —0.2) 


V.  Conclusion 

The  method  proposed  here  performs  the  analysis 
of  scattering  problems  by  the  use  of  a MoM  tech- 
nique numerically  implemented  by  wavelet  expan- 
sion. In  particular  the  use  of  the  integral  operator 
in  the  wavelet  domain  makes  it  possible  to  obtain 
the  unknown  current  without  the  use  of  quadra- 
ture formulae,  and  the  high  interpolating  proper- 
ties of  the  chosen  wavelet  basis  gives  results  with 
low  CPU  time.  The  good  quality  of  the  results  is 
demost rated  by  comparisons  with  standard  MoM 
and  FEM  computations. 
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